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We show variationally that instability of the ferromagnetic state in the Hubbard model is largely 
controlled by softening of a long-wavelength spin-wave excitation, except in the over-doped strong- 
coupling region where the individual-particle excitation becomes unstable first. A similar conclusion 
is drawn also for the double exchange ferromagnet. Generally the spin- wave instability may be 
regarded as a precursor of the metal-insulator transition. 
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I. INTRODUCTION 



Recently the electron correlation effect in the strong coupling Hubbard model around half filling has been intensively 
^"5 ' investigated. Since Nagaaka showed the existence of itinerant ferromagnetic ground state in the limit U ^ oo and 
^\ 72 ^ la several authorsoQ attempted variational estimate of the stability of the ferromagnetic state in this limit. On 
the other side, there are works trying to investigate the metal-insulator transition from the ferromagnetic side. In 
particular, the two-body problem of a particle^hole pair in the half filled band was treated as an exactly solvable case 
' of the Mott transition by several authors .ETllil To describe physically relevant situations in this approach, one must 
, treat the many-body problem of the particle-hole bound states ..Fbr .example, one may use the BCS-type mean-field 
■ approximation when the ground state is magnetically ordered.EilE3il3 

In this paper we discuss stability of the Nagaoka ferromagnetic state, with the problem of the metal-insulator 
f — transition in mind. We show that the ferromagnetic state in the over-doped region of the strong coupling Hubbard 
model is destabilized by the individual-particle excitation, as Shastry et al. noted.cl However it is fou nd that in almost 
all the other region the instability is controlled by softening of the spin-wave stiffness. In |III A| , we give a phase 
diagram showing this feature on the basis of a variational trial state. We estimate an upper bound for k, which is 
defined by k = zt/UcX for the ferromagnetic threshold in the limit [/ — > oo and x — I — n 0. Stability of the 



double exchange ferromagnet is discussed in III B| . Incidentally, in Appendix we discuss that the antiferromagnetic 



Heisenberg model can be reproduced by using the results of the two-body problem, or from the spin wave in the 
insulating ferromagnetic state. Mathematical details are given in Appendix [~ 

o 

O ■ II. VARIATIONAL DESCRIPTION OF THE SPIN- WAVE EXCITED STATE 



A. Random phase approximation 



First investigate the ferromagnetic state. 



iF)^ n 4tIo), (1) 



of the Hubbard model, 

i,j,cj i 

= Y£kcl^Cka + uYni^n,i. (2) 

k(T i 

By way of illustration, here we present results obtained in the random phase approximation (RPA) first. To this end, 
we may use the following trial state for the spin-wave excited state; 

\^,)-bl\F), (3) 

where 
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and L denotes the total number of lattice sites. Then we obtain 

(F|fejH,fet]|F) = i^n,(£fe+,-£, + C/,,)|/,(fc)p 
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and carrier density n is defined by 



Taking functional derivative of 
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with respect to fq{k), we obtain 



^nkfq{k). 



Summing nkfq{k) over k, we have the eigenequation 
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Substituting fq{k') 



Jkk' 



in Eqs. 



Ek + Un- 
(0) and (||), we obtain 

'7q(fc) = £k+q - £k + Un, 
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for the energy of the particle-hole continuum. The bound state (spin wave) energy ujq is given as a solution of Eq. (|T^) . 
In particular, for a tight-binding dispersion in a square lattice, the results are shown in Fig. |^ for H = {q, q) {0 < q < tt) 
for various values of n. The spin-wave part of Fig. |l| is shown in Fig. ||. The result for n = 1 reproduces the two-body 
result given in Appendix Discussion based on Eq. ( pT| ) is equivalent to the random phase approximation which 
properly takes into account the two-body correlation effect of the particle-hole ladder. In all of the cases shown in 
the figures, the ferromagnetic state is unstable to the spin-wave excitation. 

From these results, we observe several points, (i) ujq and riq{k) are separated by energy gap of order U. (Fig. |^.) 
(ii) The band width of ujq becomes narrow and (iii) the minimum of ujq moves away from q = Q = (tt, tt) as density 
of holes increases (Fig. ^). For example, the minimum of ujq for U/At = 4 and n = 0.9 lies at q = (tTjO.Stt). With 
respect to (iii), we are led to the following speculation; beyond spin- wave instability, the spin wave with momentum 
q = Qmin which gives the minimum Uq < will be set to populate the ferromagnetic state, resulting in Bose-Einstein 
condensation of the boson b^q^^^- In particular, around half filling n < 1, the resulting state will be the commensurate 
Neel ordered phase with q = Q. Upon doping n < 0.93, e.g. for U/At = A, the resulting phase will become the WHrjal 
state with incommensurate modulation q 7^ Q. Qualitatively this is consistent with the result of recent studies. EJ~t£l 
Moreover, further hok-daping {n < 0.82 for U/At — 4) stabilizes the ferromagnetic state, just as concluded from the 
mean-field treatment .Ej'Ej However, the latter point on the stability of the ferromagnetic state as well as the above 
point (i) are shown to be modified by improving the approximation. 
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B. Improved trial state 



Next we consider the following creation operator to improve the trial state created by Eq. (|^); 

bl = ^Yl - ''Ml (sin^ + cos9c,^cl^) c.^e^'^^^. (13) 

The wavefunction of this form was first used by Roth@ and also adopted later by Shastry et al.i to investigate stability 
of the Nagaoka ferromagnetic state. However, since the spin wave spectrum for general q and finite U derived from 
(|l3|), which turns out to be important for our purpose, has not yet been thoroughly investigated, we derive results by 
ourselves from the outset. Mathematical details are deferred to Appendix ^ Below we show only results to compare 
them with those given in the last subsection. 

For the square lattice, the bottom of the continuum ?/qniin a-nd LUq for ci — {q,q) (0 < q < vr) and U/4:t — 5 are 
shown in Fig. ^, and ojq for various values of n are shown in Fig. ^. These are to be compared with Fig. ^ and Fig. H, 
respectively. As for Uq, in the slightly doped region 1 — n <C 1, the results are not modified considerably from those of 
the RPA. On the other side, the individual particle-hole spectrum presents a striking contrast as is clear from Fig. |l| 
and Fig. ^. In the improved estimate, the continuum no more has energy of order U, but it forms a flat band lying in 
the low energy region. This is because of the fact that vacancy made through the hole doping enables particles to hop 
around, though it is a quite restricted motiona As the density of hole 1 — n increases, the band width of T]q broadens 
and finally we have a vanishing binding energy Aq = Ty^min — for q which gives the lowest w^, as shown in Fig. |^. 
In the over-doped region, therefore, the spin wave for q ~ qmin cannot be regarded as a well-defined bound state. 



III. INSTABILITY OF THE FERROMAGNETIC STATE 



A. Hubbard model 

In our previous papers,0lli we discussed stability of itinerant ferromagnets on the basis of a model comprising 
degenerate orbitals. In the strong coupling limit, we observed that the instability condition derived from the individual 
particle excitation is more stringent than the condition concluded for the spin wave instability. Therefore, putting 
more emphasis on the study of the individual particle excitation than for the spin wave spectrum, we calculated the 
critical interaction [Z^below which excitation energy of the individual particle- hole pair becomes negative as a function 
of carrier density n.EZi The result for the Hubbard model was that Uc approaches a finite value as n ^ 1 both for a 
square and a simple cubic lattice. In other words, we could not prove instability of the Nagaoka ferromagnetic state 
in the under-doped strong-coupling region, even though more elaborated trial state than that derived from Eq. ( [T^ ) 
is used to estimate energy of the individual-particle excitation. On the contrary, for the Hubbard model, a simple 
argument can be given, indicating that Uc should become infinity as n ^ 1. To show this, one may consider the case 
where holes of concentration x — 1 ^ n are doped into the half-filled Hubbard model. If x is small enough, energy 
of the complete ferromagnetic state is given by —xzt per site. On the other hand, energy of the antiferromagnetic 
configuration is —2zt^/U -\- 0{x) per site. Therefore, equating these two, the critical boundary is expected to take 
the form zt/Uc x, as was noted by NagaokaJd This argument however does not tell us whether the instability is 
brought about locally (continuously) or globally. We show that it is in fact given as a local instability by investigating 
the spin wave excitation for a strong but finite interaction energy. It is noted here that a quantitative aspect of the 
above fact was addressed in a recent work by Hanisch et al,Q where emphasis is put on how far the stable region of 
the Nagaoka state can be reduced, rather than its physical origin of our concern. 

We display the insicrse of the critical coupling as a function of n in Fig.^ and Fig.|^ for a square and a simple cubic 
lattice, respectivelyO In the figures, threshold curves determined using Eq. ( |l3|) are shown. These are calculated by 
the conditions rj^in = (dashed), ujq = (long-dashed) and D — (solid line), where ?7min is the minimum value 
of r]q{k). In the region above the curves, the ferromagnetic state is absolutely unstable. Mathematical details for rjq, 
LOq apjd D are given in Appendix The stiffness co nsta nt D for the case zt/U = was shown as a function of n in 
Ref.O. In FigjHI threshold by I?rpa — using Eq. ( B24 ) is also shown. It is clear how the trial state (|l^) improves 



the result of the RPA; in the latter we cannot prove instability for any n in the strong coupling region. 

The figures show that, in the strong coupling limit zt/U = 0, the individual particle excitation (dashed curve) brings 
about instability prior to the softening of the stiffness constant (solid curve), just in accordance with our previous 
results. However the phase boundary in the other region is primarily determined by the spin-wave instability, D < 0. 
In particular, in the region U oo and — > 1, the phase boundary is of the form zt/Uc — i^x and is determined by 
the spin- wave instability. In this limit, we cannot distinguish the two results by loq — and _D = 0, which in turn do 
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not differ appreciably from -Drpa = 0. Physically we may say that the spin-wave instability determines threshold for 
ferromagnetism in the region where the metal-insulator transition is likely to occur. A similar behavior is observed 
also for the double exchange ferromagnet as shown in the next subsection. 

As for K = zt/UcX, our calculation gives k, — 1.08 for the square lattice [z — A). Fai- the simple cubic lattice {z — &), 
we obtained k = 2.08, which is better than n = 3.96 of Richmond and Rickayzen,EI who estimated k by assuming a 
flipped spin to stayr^at a single site, i.e., not to hop around in a lattice as the spin wave does. For reference, we cite 
Nagaoka's estimatetl k — 1.47 for a simple cubic lattice. This value, however, is not to be compared with our result 
since the former does not have a variational significance. Our result, being based on the variational treatment, sets 
the exact upper bound for the true value of k. 



B. Double exchange model 



Next, we consider the ferromagnetic Kondo lattice (double exchange) model, 

where we assume a positive coupling J' > 0. This model with Sf = 3/2 is often used to describe the lanthanum 
manganese oxides Lai-jjA^-MnOa, where A is a divalent ion such as Sr, Pb or Ca. To create a trial state for the spin 
wave excitation, we use the operator 

^ 7f E ^i^"-^ - ((^i^^T + SJ,) sin e + (4^c.T + S-^)c\^c,^ cos o) ^'^^^ , (15) 



which reduces to Eq. ([l3| ) when Sf — 0. The expression for d — Q, when U = J' — 00, was previously treated by us.0 
As an example, we show threshold for the instability of the double exchange ferromagnet, the ferromagnetic ground 
state of (14). In Fig.^, we show 6tn/gc determined by = as a function of carrier density n (solid curve) where 
g = JSf + Un represents the mean- field exchange splitting of the modeLiH). We assumed a tight-binding band 
in a simple cubic lattice. In the figure, we juxtaposed our previous resullO (dashed curve), which was obtained by 
investigating the instability of the individual particle excitation; shown as the dashed curve in the right figure of Fig. 9 
of Ref.im. Experimentally, the itinerant ferromagnetic state is observed only in a restricted range 0.2 ^ x = l — n ^ 0.5. 
On the assumption that the model ( |l4|) adequately describes the manganese oxides, we may conclude that the observed 
result is explained by assuming U /6t = g/Qtn ~ 5, i.e., tJ /W ~ 2.5 with the band width W = 12t. In particular, our 
results indicate the spin wave instability in the under-doped regime, while in the over-doped regime the instability is 
controlled by the individual particle excitation with wavenumber q ^ kf. Generally the spin wave tOq is made unstable 
first at the momentum q ^ 0-,= (tt, tt, tt). Therefore the resulting phase beyond the instability is expected to be the 
incommensurate spiral stateJ23 as in the case of the Hubbard model. 

In FigJ|, we show the Sf dependence of the phase boundary Atn/g in a square lattice. Solid curves are determined 
hy D = 0. As above, for the instability due to the individual-particle excitation (denoted by long-dashed curvesi, 
we adopted the more stringent condition ek=oi = ^/ the basis of our previous results (7-9) and (7T0) of Ref.EZi, 



than that concluded from the expression derived from ([15|), the counterpart of Eq. (B6). The figure shows that the 
threshold near n ^ 1 is determined by D = 0, as in the case of the Hubbard model. The D = portion of the phase 
boundary increases as a function of Sf. In particular for Sf = 00 in the square lattice, the boundary is exclusively 
determined by the condition D — down to n = 0, while we found that the solid curve D = and the dashed curve 
cross each other around {n,6tn/g) = (0.5,0.5) for 5*/ = 00 in the simple cubic lattice (not shown in Fig.|^). Note 
that in the case Sf = 00 t he bo undar y is th e same as the result of the random phase approximation. (See I?rpa = 
of Fig.^, as well as Eqs. ( B24 ) and ( B25 ).) We see that all the boundaries in Fig.^ approach to a finite value as 



n ^ 0. This is a specific feature of two-dimensional lattices. Generally for Sf ^ 0, the parmeter 9 in ( |15[ ) determined 
variationally increases from zero to 7r/2 as n decreases from 1 to 0. Thus one can show that the condition D = gives 



Atn/g — I/tt in the limit n ^ for the square lattice, using Eqs. (B24) and (B25). On the other hand, in the simple 
cubic lattice, we have 1/U — + in this limit (Fig.^. Finally we note that the physically relevant situation Sf — 3/2 
lies just in between the classical {Sf = 00) and the quantum {Sf — 1/2) limit. 
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IV. DISCUSSION 



In the previous section, we showed that the ferromagnetic state in the Hubbard model as weU as in the double 
exchange model is unstable to the spin-wave excitation in the under-doped region, while it is unstable to the individual- 
particle excitation in the strong-coupling over-doped region. Generally the spin wave mode cjq may take a mimimum 
at finite wave vector q = qo- This minimum is interpreted to indicate the potential spiral-spin correlation contained in 
the ferromagnetic configuration, which becomes conspicuous as the excitation gap Wg,, approaches zero. We found that 
the softening of the long- wavelength spin wave mode D closely follows the gap collapse ujqg — > 0. Thus we may 
conclude that the spin-wave instability is a precursor of the metal-insulator transition. It is physically plausible that 
instability of the ferromagnetic state in the region where the metal-insulator transition is likely to occur is dictated by 
the spin wave instability, since the spin wave is nothing but a particle-hole bound state in the ferromagnetic vacuum 
and, on the other side, we may interpret the Mott insulator as composed of the particle-hole bound states (Appendix 
^). Therefore we expect that the above conclusion is generally valid; for example, instability of ferromagnetism in the 
situation appropriate to metallic nickel will be caused by the individual-particle excitation since the filling n ~ 0.2 per 
band is far from being critical for the metal-insulator transition. In this respect, we note that the Hubbard model as a 
model for an itinerant ferromagnet is a rather exceptional case, because the ferromagnetic state if any can be realized 
only around half-filling (Nagaoka limit). This is the reason why we found a dominant role played by the spin- wave 
instability in Fig.|| and Fig|6[ 

Finally let us speculate on the approach to the metal-insulator transition from the ferromagnetic side. We could 
derive the antiferromagnetic Heisenberg model from the spin-wave dispersion of the ferromagnetic insulating state, 
and observed that, unlike the individual-particle excitation spectrum, the spin-wave dispersion itself does not change 
drastically upon hole doping (Fig. |l| and Fig. . Then if the spin wave as a particle-hole bound state is robust even 
when slight amount of holes are doped, the paramagnetic phase realized after the spin-wave instability may show 
anomalous metallic properties. This kind of consideration make sense just around half filling only where the lowest 
energy spin- wave state can be regarded as a well-defined bound state (Fig. |^). Although it is an interesting problem 
to consider the interacting ferromagnetic spin wave as an elementary constituent, further investigation on this point 
requires complicated calculation which is far beyond the scope of this article. 
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APPENDIX A: TWO-BODY PROBLEM AND THE HEISENBERG MODEL 



Let us consider the simplest case of a particle-hole pair in the ferromagnetic band insulator, 

|F)^n4TlO>- (Al) 

k 

Then the eigenequation for ujq, which is the eigenvalue of the state (^, is given by 

Ek + q-Ek + U -UJq U' 

This is obtained as a limit n 1 of Eq. ([l|) , and is the exact result of the two-body problem. It is also derived 
as the condition for the particle-hole ladder to have a pole. We note that the bound state solution LOq is nothing but 
a spin wave in the ferromagnetic vacuum ([A1|). On the other side, the internal structure of the bound state fq{k) is 
given by 

fq{k) = 4f ■ z TTTT ' (^3) 

VL £k+q ~ Sk + U - UJq 

without loss of generality. 

We investigate a general case of the tight-binding band, for which Sk is given by 
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AkS 



(A4) 



where the sum is taken over nearest-neighbor vectors S. To the accuracy of order 0{t/U), from Eq. (A2) we obtain 
an expression for Wg, 



k 



where Eq. ( |A4| ) is substituted. The case n = 1 of Fig. || and Fig. ^ can be wcU fitted by this expression. 
As for fq{k), we have 



U 



to the accuracy of order 0{t/U). We then obtain 



(A5) 
(A6) 



(A7) 



(A8) 



and 



(A9) 



The result Eq. ( A9) shows that the boson b\ for non-zero t/U has an internal structure extending to neighboring sites 
of the site where it is created. Physically this structure is interpreted as a singlet cloud formed with the neighboring 
sites. 

Energy of the localized boson is given by 



{¥\h[H,b\]\F) _ 2zt^ 



(F|6.6j|F 



U 



(AlO) 



where z is a coordination number, z = X^jl- The result Eq. (AlC) is obtained also as the center of gravity of the 
band LUq, Eq. ( A6), 



2zf_ 



(All) 



From the form of the structure, Eq. (A£), Eq. (AlC) is interpreted as — z • 2t'^ /U, i.e., as a sum of energy of singlets 
formed with the z nearest neighbors. The factor 2 is due to the two processes due to a particle and a hole hopping, 
which are represented in the two terms in the parenthesis of Eq. ( A9). 

Now let us introduce the creation and annihilation operator of the hard-core boson bl and bi in place of bl and 6^; 



bl^bl. 



(A12) 



The bound state created by b} extends only to the nearest-neighboring sites of i, as indicated from Eq. (A9). As a 
result, total energy of two localized bosons differs from 2L0i only when they are in the nearest-neighboring sites, when 
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the total energy amounts to —4:{z — l)t^ /U . Increase by an amount 2 J = At'^/U from 2uji is due to overlap of the 
singlet cloud of the two bosons. In the hard-core boson picture this must be regarded as an interaction energy, i.e., 
the interaction part of the Hamiltonian for the hard-core boson is given by 



v^2jj2bikbyb,, 



(A13) 



where the sum is taken over the nearest-neighbor pairs. On the other side for the hopping part, the one-body energy 
ujq is given by Eq. (AC). In terms of J = 2t^ /U we rewrite it as 



1 



q s 



(A14) 



Instead of the hard-core boson, we can equivalently use the quantum operator for the spin S — 1/2, which are 
defined by 



S^^ = ]^- ilk, 

s+ s- = bl 



(A15) 



Then as the effective model to describe the half-filled Hubbard model in the strong cou pling regime, we can reproduce 
the antiferromagnetic Heisenberg model in terms of these spin operators: Putting Eq. ( A13 ) and Eq. ( A14 ) together, 
we obtain 

H ^f + V 



2jEfs,.S,-i 



(A16) 



Similarly as above, as an effective model for the case when holes are doped, the following Hamiltonian is suggested; 



^,3 



(A17) 



The first term describes the hopping process of a doped hole, fj = Ci|(l — blbi), and the second term takes into 
account the hopping of the boson when its neighboring sites are vacant. The Hilbert space of Eq. (A17) is spanned 
at each site by | ti), I li) = bl\ Ti) a-nd vacancy |0i) = Ci|| ti)- The fermion operator Cii must thus be operated on 
the sites where there is no boson. To take this into account, one may use // instead of c^^. 



APPENDIX B: CALCULATION FOR THE IMPROVED VARIATIONAL STATE 



We calculate uj„ for the creation operator. 



(Bl) 



2J 
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To obtain results, we may follow the same procedure as given in II A. First we obtain 



k k k k,p 



J E ^k{ek+q - £k)fk ) ( E 



L 

Jk ~^ Jk 
fk fk 



n - n 

fk fk 



(B2) 



and 



/fe -> 4* 

fk fk 



+ — ^ 2_^nkfkfk 



fk -> 4* 
/fe ^ /fe 



(B3) 



Here we denoted fq{k) simply as fk- Bracketed expressions me an to repeat the preceding terms with the replaced 
fimctions as indicated in the brackets. A^, and Tkpq in Eg. (B2) is the case 5/ = of the expression defined in (4-6) 
and (4-7) of Ref.EZl. In the tight-binding dispersion Eq. (A4), the former is given by 



Afcg-|eg|-(l-n)efe+ (1-n) 



£k- 



-9' 



where 



= Y E "^fe^fe- 



(B4) 



(B5) 



For the wavefunction Eq. (|13|), one may replace fk and fk in the above expressions by fk sinO and ffc cos 6, respec- 
tively. A function fk and a parameter d have to be fixed so as to minimize [tiM . We note that Eq. (|l^ ) for 9 — 77/2 
gives Eq. (Q). Moreover, Eq. ( p^ becomes the trial state for the case U = ooHu by assuming 9 — 0, since in this case 
the variational state does not depend on the interaction energy U at all. For a finite value of U, the parameter 9 takes 
a value in the range < 9 < tt/2. 



1. Individual particle excitation 

To derive excitation energy of the individual particle-hole pair, we may set fk' — 6kk' in Eqs. 
calculate 

^ _ (F|6Jg,bt]|F) 
(F|6,4|F). 

Akq cos^ 9 + {sk+q — £fc + Un) sin^ 9 + 2(1 — n){ek+q — £fc) sin 6* cos 6* 



Then we obtain 



77,(fc,0) 



and 



to 



(1 - n) cos2 6* 2(1 - n) sin cos 6* -h sin^ 9 
For r\q{k) this expression should be minimized with respect to 6*, and ?7min is defined as a minimum of 'r]q{k). 



(B6) 



2. Spin wave dispersion 

To min imize Uq with respect to fk, we take functional derivative dujq/df* after replacing fk and fk in Eqs. 
and (B3) by fk s\n9 and fk cos 9. Then we obtain an equation. 
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^Akqfk + ^kpqff}j COS^ 9 + (^{Sk+q - £fc + Un)fk ^ ^''f'^ ^"^^ ^ + 

^/ £fc+^_£fc ^^^^^ + 2^^k{£k+q - £k)fk + 2(1 - n){ek+q - £k)fk \ sin 6* cos 
^ k k ' 

= ujq ^(1 - + i ^ nfe/fe j cos^ 6 + fu sin^ 6* + X! ^''^'^ + ^ '^^-^'^^^ ^ '^"^ 

We investigate the tight-binding model in a square [d = 2) and simple cubic lattice (c? = 3), for which 

1 

Efc = --^cos(/ci). 

2—1 



(B7) 



(B8) 



Here and below we set zt — \ where z ~ 2d. Furthermore, we are interested in the dispersion LOq for q along the 
diagonal of the Brillouin zone, i.e., for q — and q = {q,q,q) (0 < q < n) in the square and the simple cubic 

lattice, respectively. Then we can cast Eq. (B7) into the following form. 



(B9) 



where 



D[k,q) = (Afcg - (1 - ii)LOq)cos^ + {sk+q - £k + Un- Wg) sin^ 6* + 2(1 - ii){£k+q - £k- ujq) sin6'cos6', (BIO) 



Piil) = 7 X! ^kiSil - £kSi2 + £k5a)fq{k), 



(Bll) 



Ni{k,q) = (ujq - \eg\£q - (1 - n)ek+q) cos^ + U sir? + {2ujq - Sk+q + £k)sm9cos0, 



N2{k,q) = - {\eg\ek+q + (1 - n)eq) cos^ 6* - (1 + Sq) sin ^ cos ^ 



Nsik^q) = - {\eg\£k+q + (1 - n)eq) 



cos t) — v„ sm U cos ( 



Here we introduced 



1 

Efc = - ^ sin(A:i), Vq = sm{q). 



Solving Eq. (B9) for fq{k) and substituting the result into Eq. ( Bll ), we get the eigenequation 

det {A,j{q)-6^j) = 0, 

where 



N,{k,q) 
D{k,q) ■ 



(B12) 
(B13) 
(B14) 

(B15) 

(B16) 
(B17) 



To obtain the spin wave energy ujq, we must minimize the solution ujq{0) of Eq. (B16) with respect to 6. 



3. Spin wave stiffness constant 



The spin-wave stiffness constant D is defined by 



(B18) 



In the long wavelength limit g — > 0, we can expand fq{k) with respect to q and derive D analytically from Eq. (|B2| ) 
and Eq. (B3). Below we give only lie resulting expressions. 

As shown in our previous paper ,113 the results for the Hubbard model are obtained as a special case of those for the 
Kondo lattice model (p^). Using Eq. (p^), we calculated D for this general model; 



D{e) = Do- SD{e), 

where Dq = \eg\/2Sz in terms of twice the spontaneous magnetization 25* = 25*^ + n, and 



dD{d) 



I ((1 - n) cos2 e + sin^ 9 + {2 - n) sin 9 cos 9)' 



25 ( 1 + I eg I z/ cos2 6'/2) ( 1 + 2 sin 61 cos ( 
In these expressions, /, and Ako are given by 



(> 0). 



L ^ Ako cos^ 9 + gs'm 9 



(B19) 



(B20) 



(B21) 



dkx 



and 



Afeo = {2Sf + 1) (|e<,| - (1 - n)ek) + ((1 - nf - le^p) (> 0). 



(B22) 



(B23) 



A parameter g in Eq. ( B21 ) is defined by g — Un + J'Sf, and represents the mean-field exchange splitting of the model 
(|l|). We must minimize D{6) to obtain D. To reproduce the result for the Hubbard model, we may set Sf =jQ— .The 
particular case of the above results, i.e., in the strong coupling limit g ~ oo (9 = 0), was obtained previously.uO For 
the Hubbard model, the result of the random phase approximation -Drpa is obtained by setting S f ~ and 9 = tt/2, 



-Drpa 



zn 



1 1 



Un^ L 



(B24) 



A similar expression is obtained in the limit S —> oo where we should assume 9 = t:/2 to keep / and 2S5D{9) finite. 
Then we obtain 




1 1 



E 



nkvl^ 



(B25) 



In the main text, we regarded the approximation scheme giving the result corresponding to (B24) and ( B25| ) as the 
random phase approximation. [EJie former result (B24) is well known for the Hubbard model. The latter was recently 
obtained by the other method.Eil 
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FIG. 1. Dispersion of the spin wave (ujq/it) and the bottom of continuum (»7gmin/4t) along {q,q) (0 < g < tt) in the RPA. 



FIG. 2. Dispersion uj{q) for U/U = 4 in the RPA. 



FIG. 3. Dispersion of the spin wave {ujq/4:t) and the bottom of continuum (rjqmin/it) along (g, g) {0 < q < n) calculated 
with the improved trial state. 



FIG. 4. Dispersion uj{q) for U/4:t = 4 for the improved trial state. 



FIG. 5. Threshold for the stability of the ferromagnetic state in a square lattice. 



FIG. 6. Threshold for the stability of the ferromagnetic state in a simple cubic lattice. 



FIG. 7. Inverse of the critical coupling Gtn/gc as a function of n for the Sf = 3/2 ferromagnetic Kondo lattice model in a 
simple cubic lattice. The ferromagnetic state is unstable outside the region denoted by 'F'. 



FIG. 8. 4tn/(?c as a function of n for the ferromagnetic Kondo lattice model in a square lattice for Sf = 1/2, 3/2 and oo. 
Solid curves are determined by D = 0. Instability of the individual-particle excitation occurs above the long-dashed curves. 
The ferromagnetic state is unstable outside the region denoted by 'Ferromagnetic'. 
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